arXiv:1501.00757v2 [hep-th] 17 May 2015 


Prepared for submission to JHEP 


Renyi entropy of locally excited states with thermal 
and boundary effect in 2D CFTs 


Wu-Zhong Guo, a b Song He, 6a 

a State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese Academy 
of Science, Beijing 100190, P. R. China 

b Yukawa Institute for Theoretical Physics, Kyoto University, Kitashirakawa Oiwakecho, Sakyo-ku, 
Kyoto 606-8502, Japan 

E-mail: wuzhong@itp.ac.cn, hesongl7@gmail.com 

Abstract: We study Renyi entropy of locally excited states with considering the ther¬ 
mal and boundary effects respectively in two dimensional conformal field theories (CFTs). 
Firstly, we consider locally excited states obtained by acting primary operators on a ther¬ 
mal state in low temperature limit. The Renyi entropy is summation of contribution from 
thermal effect and local excitation. Secondly, we mainly study the Renyi entropy of locally 
excited states in 2D CFT with a boundary. We show that the time evolution of Renyi 
entropy is affected by the boundary, but does not depend on the boundary condition. 
Moreover, we show that the maximal value of Renyi entropy always coincides with the log 
of quantum dimension of the primary operator. In terms of quasi-particle interpretation, 
the boundary behaves as an infinite potential barrier which reflects any energy moving 
towards it. 
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1 Introduction 

Many kinds of observables can be defined in Quantum field theories (QFTs). When we 
study global or non-local structures, entanglement entropy (EE) or the entanglement Renyi 
entropy (RE) are very helpful quantities. For a subsystem A, both of them are defined 
as a function of the reduced density matrix pa■ The reduced density matrix pa can be 
defined from the original density matrix p by tracing out the subsystem B which is the 
complementary of A. 

One may be curious about whether there is a kind of topological contribution in entan¬ 
glement entropy even for gapless theories, particularly for conformal field theories (CFTs). 
For example, topological properties can be qualified by computing topological contribu¬ 
tions in entanglement entropy called topological entanglement entropy [1]. In this paper, 
we focus on extracting such kind of topological quantity from both Renyi entropy and von- 
Neumann entropy of locally excited states in two dimensional rational CFTs with thermal 
effect and boundary effect. 

The n-th Renyi entanglement entropy is defined by S A ^ = log Ti'[/9^]/(l — n) formally. 
The limit n —>• 1 coincides with the von-Neumann entropy. This is standard replica trick 
method to calculate the entanglement entropy. The difference of S A ' between the locally 
excited states and the ground states with introducing thermal and boundary effects are 
main interest in this paper. The difference is denoted by A S A \ Replica approach calcula¬ 
tions of A S A for states excited by operators have been given in [2, 3, 5]. We will closely 
follow the construction in [3][4][5]with introducing thermal and boundary effect. 
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We would like to review the thermal effect and boundary effect respectively. There 
are many studies about thermal effect in 2D CFTs. For thermal states, EE is not a good 
measurable quantity, which is contaminated by the thermal entropy of the subregion. In 
high temperature limit, the EE will be dominated by thermal entropy. To reveal the 
quantum entanglement of system with thermal effect, one should identify the thermal 
contribution and other contributions of EE. Ref. [6] conjectured universal form of correction 
of EE in any quantum system with mass gap. Ref. [7] provided the form of the coefficient 
of such correction in 2D CFT. To generalize studies in [7] to higher dimensions, Ref. [8] 
and Ref. [9] considered thermal corrections to the entanglement entropy on spheres. On 
the other hand, the dynamics in 2D CFTs with a boundary have many new features 
comparing with 2D CFTs in full complex plane. The original works have been done by 
Cardy, who discussed surface critical behavior of correlation functions [10]. Ref.fll] studied 
the constraints on the operator content by imposing by boundary conditions and also the 
classification of boundary states in terms of the modular transformation. In [12] [13], the 
concept of boundary operators have been introduced. Ref. [14] showed that the resulting 
set of boundary conditions to be complete. There are also nice correspondence called as 
AdS/BCFT proposed by [15] [16]. The boundary effect can be also studied holographically, 
which is beyond the issue considered in this paper. 

In 2D rational CFTs, the authors in [4] get an amazing result for the locally excited 
states, which relates the Renyi entropy to the quantum dimension of the primary operator 
which is kind of topological quantity. In this paper, we generalize the previous study [4] 
on Renyi entropy with thermal and boundary effects. Firstly, there is a simple sum rule 
between the thermal correction and local excitation in low temperature limit. That is to 
say the total Renyi entropy are summing over Renyi entropy of local excitation and the one 
of thermal excitation in low temperature limit. Such kind of relation is similar to the sum 
rule related to the Renyi entropy in [5]. One can generalize the result to local excitation 
in pure state in 2D CFT. We make use of a different approach [6] to obtain the thermal 
correction to Renyi entropy which can be reduced to [17]. Secondly, we investigate the the 
Renyi entropy for states excited by local primary operators in the rational CFTs with a 
boundary. These boundaries introduced here do not break the conformal symmetry. Such 
theories are called as BCFTs. We show the time evolution of the Renyi entropy in 2D 
free held theory and Ising model. Then we generalize to rational CFTs with a boundary. 
The boundary changes the time evolution of the Renyi entropy, but does not change the 
maximal value of the Renyi entropy. All these cases studied in this paper show that the 
Renyi entropy does not depend on the choice of boundary conditions. In 2D rational 
CFTs with a boundary, we also show that the maximal value of the Renyi entropy always 
coincides with the log of quantum dimension of the primary operator during some periods 
of the evolution. We give the physical understanding of boundary effect, which support 
the quasi-particles explanation of the local excitation. The boundary behaves as an infinite 
potential barrier which reflects the quasi-particle moving towards it. 

The layout of this paper is as follows. In section 2, we study the thermal effect on 
the Renyi entropy of the local excited state in low temperature limit. In section 3, we set 
up the local excitation in 2D CFT with a boundary and obtain the Renyi entropy of a 
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subsystem with time evolution. We study the 2D free scalar, and Ising model as examples, 
then generalize the result to 2D rational CFT. In section 4, we devote to the conclusion 
and physical interpretations of such kinds of effects shown in this paper. 


2 Local excitation in non-vacuum states 


In this section, we would like to study the local excitation of thermal state. We consider a 

system with temperature T = 1//3 and assume the excitation is local at x = — L by primary 

operator O shown in fig. [1] . In this section, we just only consider the low temperature case 
with large f3 [8]. The subsystem A is — l < x < 0. The density matrix p(t) is 

p(t) = N(t)(0(u 2 ,0 2 )(^2\n) (n|e^ En_2 ^ e£n )O t (o;i,cDi)^, (2.1) 

where we have considered the real time evolution, e is the ultraviolet regularization, N(t) 
is fixed by normalization condition trp(t) = 1, E n are the energy of the excited states. The 
complex coordinates in u plane are listed as follows. 

uj\ = i(e — it) — L, u )2 = —i(e + it) — L, 

u>i = —z(e — it) — L, ui 2 = *(e + it) — L. (2-2) 

For later convenience, we define 


i. 


t 


x 



A 



Figure 1. This figure is to show our setup in two dimensional complex plane ui = x + it. The 
system will be triggered at x = —L and there are left- and right-moving quasi-particle at t = 0. 
This setup is same as the one in 2D CFTs without boundary [4]. 


Po(t) = tr B {0(oj 2 ,u 2 ) |0) (0| O f (wi,wi)), 

pi(t) = tr B {e~^ El ~ 20eEl O{uj2,uJ2) |1) (1| O t (wi,wi)), (2.3) 


which can be taken as the reduced density matrix related to the vacuum and first excited 
state respectively, where we normalize the vacuum energy to be zero, B is the complemen¬ 
tary part of subsystem A. In the low temperature expansion with (3E\ <S 1 


p A (t) = tr B p(t) 


Po(t) + pi(t) + ... 

tVA(po(t) + Pl(t)) + •••' 


(2.4) 
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In terms of the definition of the Renyi entropy 


c (n) _ logtrp A (t) r 
— 


1 — n 
1 

1 — n 

1 


1 — n - 


log 

log 


tr(po(t) n 

(trp 0 (t)y 

tr(po{t) n 


■(1 + 


ntr(p 0 (t) n 1 pi{t)) ntrpi(t). 


+ 


tr(po(t) n ) 
ntr{pQ{t) n ~ l px^t)) 


trp 0 (t ) 
ntrp\(t ) 


(2.5) 


0 trpo(t)) n tr(p 0 (t) n ) trpo(t) J 

When there is no local excitation, i.e., the operator 0 = 1, the result is the same as [7]. 
The second and third terms of the last line in (2.5) involve in the coupling between the 
local excitation and thermal environment. In terms of the state operator correspondence, 
one can denote |1) = |0) [7]for the excited state with energy E\. Here we 

consider the whole system has an infrared cut-off A, and l/A -C 1. Using the path integral 
language, (2.5) is 


c( n ) _ 


1 


■log 


(0 t (wi,Wl)0(w 2 ,a;2)...0(w2n,W2n))c' Tl 


+ 


1 -n"° ((Ot^o^O^,^))^)" 

ne~ l3El (O^ (cai,thi) 0 (a; 2 ,o) 2 )... 0 (cu 2 n;tU 2 n)V , (—oo)V^(+oo))(g ri 
1 - Tl {0^{ui,Qx)0{u2,U2)...0{u2n,U2n))c n 

ne~P El (O t (cui,0i)O(a;2,02)^(-oo)^(+oo))c 1 

1-n (Ot( Wl) a)i)0(a;2,a; 2 ))c'i 

where C n is the n-sheet cylinder with circumference A and 0(wi, wi)...0(w2ti, W 2 n) are the 
operators that are inserted in the suitable place in the n-sheet cylinder. The first term 
(2.6) is given by [4] for the local excitation in vacuum. We will study the second and third 
term of (2.6) in detail. 

The following transformation [19] can map the n-sheet cylinder to a cylinder with 
circumference A, 

/ g2i7T£j/A _ j \ 1/ra 

~ _ g2«7r//A J * ' 

The points u = —oo and u = oo are mapping to Z-oo = e ^ 2ml / A and z +OQ = 1 respectively. 
For simplifying our analysis, we only consider the range |u;| <C A. Otherwise, the following 
calculation will be much more complicated. But the final statement does not change 
without this approximation. Thus (2.7) reduces to 

l/n 


Z ~ 


W 

W — l 


( 2 . 8 ) 


which is same as the one that is used in [4], The points z \,..., Z 2 n are given by 

„ _ Jli-Kk/n 

~2fc+i — e ' 


~ _ „2ink/n 

%2k+2 — e 


= _ „2i-Kk/n 

Z2k+i — e ' 


= _ J2mk/n 

?2k+2 — e ' 


ie +1 

- 

L 

s^l/n 

ie + t — 

L 


l) ’ 

+ 

•<r* 

1 

t 


L ^ 1/" 

—ie +1 

— 

L 

-l) 

—ie — 
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L \ !/”■ 

—ie — t 

— 

L 

-l) 

ie — t 

— 

L 


ie — t — 
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l) ‘ 


(2.9) 
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Int>/ + LorO<t<L 


2 iLe 

Z2k+l ~ Z 2k+2 - ~ _ L ^ f _ x _ L ^Z 2k +l, 

2 iLe 

Z “+i “ * 2 ‘+ 2 ~ ~ n(t + L)(t — l + L) Z2k+1 ' (2 ' 10) 

In L < t < L + l, 

2 iLe 

Z2k ~ Z 2t+1 = -„ (< _ i)(t _ 

2 iLe 

Z2/C+1 - ~2fc+2 ^ - n ( t + L )( t _l + L f 2k + 1 ' ( 2 - n ) 

In t>Z + L or 0<t<L the second term of (2.6) is 

ne~P El (O^(ui,ui)0(u 2l uj 2 )...0(u 2n ,u) 2n )ip(—oo)'il>(+oo)) Cri 
1 ~n (0^ui,Oi)0(uj 2 ,ui 2 )...0(uj 2n ,u> 2n )}c n 

_ ne~ l3El {0 :i (zi,zi)0(z 2 ,z 2 )...0(z 2n ,z 2n )'i/j'( z-oq)-^'( z +00 )) Cl 
1 - n (0^(zi,z 1 )0(z 2 ,z 2 )...0(z 2n ,z 2n ))c 1 

71 f / / 

= -(V’ (z-oo)i> (z+ oo))ci- (2.12) 

I — n 

In the above formula, we label ip' (z,z) as the map of the operator ip(u,u) 1 . In the second 
step due to (2.10) we have used 

( 0\z i, zi)0(z 2 ,z 2 )...0(z 2n , z 2n )ip' (z — OO (z+ OO ))ci 

~ (O f (zi, Zi ) O(z 2 , z 2 ))... (O f (z 2n -! i Z 2 n—l)0(z 2n , Z 2r P)) {ip (*-—oo)^ (^+oo))- 

(2.13) 


In L <t t T T l , the correlation function can not be factorized as (2.13) directly due 
to (2.11). Following logic in [4], one can make use of n — 1 times fusion transformation 
(zi,z 2 )(z 3 ,Z 4 )...(z 2n -i,z 2n ) ->• (z 2 ,z 3 )(z 4 ,z 5 )...(z 2n ,zi). The second terms of (2.6) is still 
given by in e —> 0 

up , , 

— (z-ooW (z +00 )) Cl . (2.14) 

Here we assume that there are no nontrival correlation between O and ip for simplifying 
analysis. The third term of (2.6) in the limit e —>• 0 is 

ne ~P E i u>i)0(u 2 , u 2 )ip{— °o)'0(+oo))ci 

1 ~n (0 ^(uji,Oi)0(uj 2 ,0 2 )) Cl 

= -T~— (V ; (-oo)V ; (+oo))ci = - 7 —— e _/3El . (2.15) 

1 — n 1 — n 

1 The operator may not be a primary operator , such as the energy-momentum tensor T. A Schwarzian 

derivative term related to energy momentum tensor operator [7] [18] will present due to conformal transfor¬ 
mation. 
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The sum of the second and third term is the same as the thermal correlation [7] for the 
short interval limit. (2.6) is the summation over the thermal correction of the vacuum state 
and local excitation in the vacuum state in low temperature limit. The Renyi entropy is 
summation of thermal effect and local excitation. This relationship can be considered as a 
sum rule for Renyi entropy of low temperature thermal state with local excitation, which is 
different with the sum rule proposed in [3]. For the pure state one could also get a similar 
sum rule. The local excitation in thermal states have been also studied in [17] for free boson 
and Ising model. With a different method [7], we can reproduce the Renyi entropy [17] for 
short interval in the low temperature with taking limit e —>• 0 for general 2D rational CFT. 


3 Local excitation in 2D CFTs with a boundary 

In this section, we would like to study the Renyi entropy of locally excited state in 2D CFT 
with a boundary which preserves conformal symmetry. The global property of the CFT 
with a boundary has been discussed in literature[19] [20] [21]. As we know, the Renyi entropy 
is sensitive to correlation function. The boundary will change the correlation functions. 
And the boundary conditions also affect the correlation functions. It is interesting to check 
what will happen to the Renyi entropy for the local excitation in 2D CFT with a boundary. 


3.1 Set-up of local exciation 

We begin with a CFT with a boundary at x = 0 and the CFT is living in the range x < 0. 

We divide the this region into two parts, one part is — l < x < 0 denoted by A and the 

other is complement to the region A denoted by B. The time t vary from —oo to +oo, and 
the Hamiltonian H is well defined as a operator to generate the time evolution. 

We assume that the local excitation of vacuum is at x = — L and consider the Renyi 
entropy of the subsystem A. The time dependent density matrix can be written as 

p(t) = NO{u 2 , 0 ) 2 ) |0) (0| O t (wi,wi), (3.1) 

where the coordinates are, 

oj\ = i(e — it) — L, U 2 = —i(e + it) — L , 

wi = — i(e — it) — L, C 02 = i(e + it) — L. (3.2) 


We still make use of replica trick to study the variation of the Renyi entropy of the 
subsystem A in this section. By definition, the variance of n-th Renyi entropy can be 
calculated as 


A^ n) 


1 

1 — n 


log(0^ i (LUi,UJi)0(u2,U2)...0(u2n,U2n))B rl 


nlog{0\uJi,Ui)0(u2 1 UJ2))B 1 , 


(3.3) 


where B n is the n-sheet Riemann surface that consists of n copies of original plane x < 0 
with gluing together along —l^x^0,t = 0. With the following conformal transformation, 


2 : 


n 


UJ H - l 
UJ — l 


(3.4) 
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t 


X 



A 


Figure 2. This figure is to show our setup in two dimensional left half plane lo = x + it with a 
boundary x = 0. The system will be triggered at x = —L and there are left- and right-moving 
quasi-particle at t = 0. 


the n-sheet Riemann surface can be mapped to a disc \z\ ^ 1 which is smooth surface. The 
boundary x=0 corresponds to \z\ = 1. Furthermore, we can map disc to the upper half 
plane (UHP) t ^ 0 with the other conformal map 


€ = -* 


.Z + l 


z-1 


(3.5) 


After two conformal maps, one can make use of well known results of 2D CFT on UHP. 
Finally, the variation of the Renyi entropy (3.3) is 



_J_1 [" TT / duJk , , cfcUfc ._fr (Qt(£i, 6)0(6; 6)---O t fen-l,C2n-l)Ofen,6n))t/gP ' 

1-77-0 k=l d ^ k ({0 ^(u>i,U!i)0(u}2,02))b 1 ) n 

(3.6) 


This is the key formula in the remainder of this paper. As a simple example, we consider 
the 2nd Renyi entropy firstly. These coordinates on UHP can be expressed by the original 
spacetime coordinate as follows. 


6 = ~i 


, (gk) 1/2 + i 
<sk) i/2 -1’ 


£3 = ~i 


-(S?l) 1/2 + i 


-(SU7) 1/2 - 1’ 


- = . (SUi) 1/2 + 1 

! (§hf) 1/2 -i' 


£3 = * 


-(§S) 1/2 + i 


-(SUi) 1/2 - 1 ' 


£2 


^4 


6 


& 



(3.7) 
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One can analyze the above formula carefully in two different time regions. When 0 < t < 
L — l or t > l + L with e —» 0, 


6 - 6 - 
6 - 6 - 
6 - 6 =* 
6 - 6 ^ 


26 e 


%/ Wl -T /\/ tai — l 
2&e 

\/ tai —~Z\/ — l 

26 e 

wi — Z \/wi -f- 1 

26e 


■\/wi — Z%/wi T 1 
and the Jacobi factor of conformal transformation, 


o?6 ^6 1 t 

dui di 02 lie 12 ’ 

^6 ^ o?6 ^ 1 ^ 

do ;3 duj^ lie 34 ’ 


^6 ^ ^6 ^ _ J_d , 
dw2 2*e 12 

dis d6 1 - 

——- 64 

dui 3 duj4 He 


When L — l < t < L + l with e —>• 0, 

6 - 6 - 
6 - 6 — 


26 e 


6 - 6 - —7^ 
6 - 6 ^ - 


y/ Ul\ -fly/ UJ 1 — ^ 
26 6 

v6i t Zy6i — ^ 

26 e 


V^l — ^6^1 T ^ 

2f 3 e 


6 Wi — ly/ U)\ + Z 

and the Jacobi factor of conformal transformation, 

dfji ^ d6 _ J_£ d6 ^ d6 ^ _ J_d 

dwi dui 4 lie 14 ’ du 1 did 2 lie 12 

— — —64 
2ie^ 3 


d& ^ d%3 ^ —63 ~ 

du)2 duj3 lie ’ du 3 did 4 


3.2 2nd Renyi entropy for free boson 

We will focus on the following local operators in the free scalar held firstly, 

Ch = e^ /2 , 0 2 = -^(e^ /2 + e-^ /2 ). 
v2 


(3.8) 

(3.9) 

(3.10) 

(3.11) 


(3.12) 

(3.13) 

(3.14) 


(3.15) 


The time evolution of Renyi entropy for such operators have already been studied in [4] 
in 2D CFT on the complex plane. There are two kinds of boundary conditions for 2D 
free scalar held theory. One is f^|s = 0 called Neumann boundary condition, the other is 
(J)\b = 0 called Dirichlet boundary condition. Since the boundary condition is homogenous, 
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it is invariant under the conformal transformation. 

The image method [22] 2 is an efficient way to obtain the correlation function on UHP 
from correlation function on the full complex plane. The two kinds of boundary conditions 
correspond to different parity transformation in image method. Due to the presence of 
boundary, there are constraints on local conformal transformation, the anti-holomorphic 
and the holomorphic sectors in correlation function are no longer independent. More pre¬ 
cisely, the correlation function on the upper half plane can be expressed by holomorphic 
part of conformal block on the full complex plane with including the ‘images’ of the holo¬ 
morphic coordinates. That is to say, 


equals to 


{(t>{zi,zi)^>{z 2 , z 2 )...cl){z n , z n )) UH p 


{4>{z\)<t>{z 1 *)4>{z2)4>{z2*)--4{z n )<t>{Zn))R2, 


(3.16) 

(3.17) 


where <f> and <fi refer to the holomorphic and anti-holomorphic part of the field (f> 3 . After a 
parity transformation the anti-holomorphic part become a holomorphic field with conformal 
dimension of the original anti-holomorphic part. In terms of image method, we should 
introduce parity transformation. For the free boson the parity transformation is 


cj)(z,z) = r)(j)(z,z), 77 = ±1, (3.18) 

r] = 1,-1 corresponds to the Neumann boundary condition and Dirichlet boundary con¬ 
dition respectively. In terms of (3.6), to obtain Renyi entropy, it is necessary to know the 
two-point and four-point correlation function on the UHP. 


3.2.1 Local excitation 0 \ 


Let’s consider the operator 0\ firstly. Using the image method, we could get the two-point 
function 


{0j(w 1 ,0 I )0 1 (w 2 ,w 2 )>B 1 = 


UHP 


ti < L 


- holo 


( 44) 1/4 (4c 2 ) 1 /*' 


(3.19) 


where £3 = £ 3*,£4 = £ 2 *,^ = ioJi, h = h = 1/8 and 0\ is the field with parity transformation. 
The subindex ‘holo’ means that we only keep the holomorphic part of the correlation and 
set the anti-holomorphic part to be constant 4 which is determined by boundary condition 


2 In conformal field theory literature, the image method is also called double trick sometimes. 

3 For most primary fields, it is impossible to divide the field <p into holomorphic and anti-holomorphic 
part. Here we only want to express that the n-point correlation functions on UHP which are dependent on 
the holomorphic conformal blocks of 2n-point correlation functions on the full complex plane. In 2D free 
field theory, the conformal blocks of the operator 0 1 is trivial. In this paper, we will also show the image 
method in 2D Ising model and other generic 2D CFTs. 

4 We can use this rule in 2D free scalar field theory, since the conformal blocks related 0 1 are actually 
trivial. We would like to appreciate communication with Cardy on this point. 
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in general. For two-point function one could normalize the field and take constant to be 1. 
The 4-point correlation function could be obtained by similar procedure. 


( 0 \ (wi, wi)Oi (w 2 , oj 2 ) 0 \ (w 3 , w 3 )Oi (w 4 , w 4 ))s 2 
^_1 a <,i a <,i 






. holo 


(3.20) 


where £5 = £ 1 , £e = £ 2 , £7 = £ 3 , £s = £ 4 , the operator 0\ = e ir] ^ on R 2 with parity 
transformation. 

In the region 0<t<L — l or t>l + L, as (3.8) shows, the correlation function (3.20) 
can be factorized as 

(Oi(^l,fl)...Oi(^4,^4)Oi 1 '(^5,^5)...di (£s,£8))p2 

OC (Oj(^i,^i)Oi(^ 2 ,6))_R2...( 01 ^ 7 ,^ 7)01 (^8,^8))i?2 (3.21) 


in e —>• 0 limit. The image method leaves us with a constant C. To fix the constant C, we 
take the limit £1 —>• £ 2 , £1 —> £21 £3 —>• £4 and £3 —>• £4 in (3.20), one could find 

(0l(6i£i)Oi(&,£2)Ol(£3,£ 3 )Oi(£4,£4))[/i7P 

— (Oi(6 5 6)Oi(6 ! 6))c/hp(C ) |(^3,^3)Oi(^4,^4))c/pp. (3.22) 


In terms of image method, we can also obtain 

(o\ (Ci, Ci)Oi (6, &)o\ (£ 3 , Cs)Oi (£ 4 , U))uhp 


(Ol(£l,£l)...Ol(£ 4 , ^4)0(6,6)01' (£ 5 , C5)---Oi(C8, C 8 ))p 2 

c 


t. 


. holo 


(C12C34C56C78) 1 / 4 ’ ^ ^ 

where C is constant. Comparing (3.22) with (3.23), we can fix C = 1 which is consistent 
with the normalization of the two point correlation function. 

For n = 2, the variation of the Renyi entropy (3.6) is 

l°s I"TT(^fe)-/ l (^fc)-/ l (^l(h ! 6)0(6,6)01(6,6)0^4, U))uhp- 
d£k d£i. 


a = -: 


k =1 


= - log 


(C 12 C 34 C 56 C 78) 2 


)l/8_ 


(( 0 t(o; 1 ,a; 1 ) 0 (a; 2 ,a; 2 ))s 1 ) 2 
(4c 2 ) 1 / 2 


= 0 , 


(3.24) 


(2e) 8 ' (C12C34C56C78) 1 / 4 ! 

where we have used (3.19), (3.20),(3.23) and Jacobi factor (3.9)(3.10). 

In the other region L — l < t < L + l, we can not factorize the correlation function 
as (3.21) directly. For 2D free scalar theory, the situation becomes much simpler. The 
correlation function could be expressed as [22] 


<e ; 


iaicf) e ia n <l>\ 


m 

Kj 




(3.25) 


- 10 - 














with the neutral condition a\ + + ...a n = 0 and z tJ = Zi — Zj,Zij = Zi — Zj. Thus 


(0}(6, ?l)---Ol(^4, C4)0(^2,6)01' (£ 5 , ?5)---Ol(^8, ?8))p2 


t, 


- holo 


1 


(^41^23^56^78) 1/74 ’ 


(3.26) 


( 2 ) 

where we have fix the constant to be 1. To get A S y A we only need to change 2 ■<->- 4 in 

(3.27) 


(3.24). Using (3.12)(3.19)(3.26) and Jacobi factor (3.13)(3.14), we get 

= 0, 


a 4 2) = -log 


^(6462^5668)^^ 1/8 _ ( 4 e 2 ) 4 / 2 


( 2 e) 


( 64 ? 32 ^ 56 ^ 78 ) 1 / 4 - 


( 2 ) 

in L — l <t < L + 1. To close this subsection, one more thing should be noted that A S\ 

does not depend on the choice of parity transformation. Actually in our above calculation 

we do not use the value of r], which is related to the boundary condition. The leading order 

( 2 ) 

of the correlation function is same. A S A is always zero for operator 0\. To see the effect 
of the boundary, we should consider more complicated example. 


3.2.2 Local excitation O 2 

t A (2) 

O 2 is a linear combination between 0± and 0[. We can calculate A S A with following the 
logic in previous section. The two-point correlation function of O 2 , 


{01(€i,£i)02(£2,Q)uhp = x (0\(£i,£i)0\(£, 2 ,Q) UH p + (Oi(^ 1 ,QOi^ 2 ,Q) UH p 


1 


it 


1 1 


+(Oi(£i, £i)0\(£,2, £2))uhp + £i)Oi(£ 2 , £ 2 )) UH p 

1 _ 1 
(4) 1/4 (4) 1/4 " WW 1, 


(3.28) 


in the limit e —>• 0. Here we have used the parity transformation (3.18) related to the 
boundary condition. We have set anti-holomorphic parts to be 1. 
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The four-point correlation function on UHP is 


(02 (£1 > £1) • ■ - 02 (£4, £4)) UHP 


(01 (£11 £1) 01 (£2, £2 ) 01 (£3, £3 ) 0\ (£4, £4 ) ) UHP 


T (Oi (£1 1 £i)Oi (£2 j £2)0} (£3? £3)0! (£4; £4)) c/pp 
+ (01 (£1; £1) 01 (£21 £2) 01 (£3, £3 ) 01 (£4, £4) } UHP 
+ ( 0 i(£ 1 )£i) 0 i(£ 2 ,£ 2 ) 0 i(£ 3 ,£ 3 ) 0 i(£ 4 ,£ 4 )}c/HP 
+( 0 i(£i)£i) 0 i(£ 2 ,£ 2 ) 0 l(£ 3 ,£ 3 ) 0 i(£ 4 ,£ 4 ))c/pp 
+ (01 (£1; £1) 01 (£21 £2) 01 (£3, £3 ) 01 (£4, £4) } UHP 
T ( 0 i (£1 1 £1 ) 0 i (£2) £2)0} (£35 £3)0! (£4^ ^4)} U H P 

T (Ol (£l) £l)Ol (£2) £2)0} (£35 £3)0! (£ 4 ) ^4:')'} U H P 

+ (01 (£1; £1) 01 (£21 £2) 01 (£3, £3 ) 01 (£4, £4) } UHP 
T ( 0 i (£11 £1 ) 0 i (£2 j £2)0! (£3? £3)0! (£4? £4)') UHP 
+ ( 0 l(£l,£l) 0 l(£ 2 , 6 ) 0 l(£ 3 ,£ 3 ) 0 i(£ 4 ,£ 4 ))pPP 
T(01 (£l j £l) 0 l (£2 j £2)0! (£35 £ 3)01 (£ 4 ) ^ 4 ))[/ffP 
T (Oi (£11 £i)Oi (£21 £2)0! (£3? £3)01 (£4; £4)) c/pp 
T (Oi (£11 £1 ) 0 j (£2 j £2)01 (£3? £3)01 (£4; £4)') uhp 
T ( 0 i (£1 1 £i) 0 i (£2 j £2)01 (£35 £3)01 (£4^ £ 4 )}c/pp 
T (Oi (£1 j £1 ) 0 i (£2 j £2)01 (£35 £3)01 (£4^ ^4)) u h p 


(3.29) 


We could use the result of four-point function of Oi which has been studied in section 
(3.2.1). From (3.18)and (3.17) one could see that the four-point function is dependent on 
the boundary conditions or parity transformations. We will calculate Awith Neumann 
and Dirichlet boundary condition respectively. 

1. For the Neumann boundary condition, i.e., 77 = 1. 

In (3.29), terms containing equal number of 0\ and 0\ will survive due to the neutrality 
condition. Thus there are 6 terms making contribution to the 4-point correlation function. 
In the region 0 < t < L — l or t > L +1, (3.29) can be expressed by factorized form on the 
R 2 . For example, the first term in (3.29) as a leading term is 


(Oi (£1, £i)Oi (£2, £2)01 (£3, £3)0! (£4, £4 ))uhp 
(Oi (£i)£i)ol (£2, £2)01 (£3, £3)0! (£4, £4)01 (£5, £5)0! (£ 6 , £ 6 ) 0 i (£7, £7)0! (£ 8 , £s)) 

(Oi (£1! £i)oI (£2, £2)) (Oi (£3, £3)0} (£4, £4)) (Oi (£5, £5)01 (£ 6 ,63)) 


. holo 


( 0 i (£7, £7)01 (£ 8 , £s)) 


t/ 


1 


1 

rsj — 

1 

e 


J holo (£l2£34£56£78) 1 0 

in the limit e —>• 0. The second terms as sub-leading term is 

Pint 


(3.30) 


(Ol (£l> £l)Ol (£2j £2)0} (£3j £3)0} (£ 4 ) £ 4 ))[/ffP = (Ol(£l,£l) 0 l(£ 2 ,£ 2 ) 
0 i (£3, £3)0! (£4, £4)01 (£5, £5)01 (£ 6 , £ 6 ) 0 i (£7, £7)0! (£ 8 , £s))p2 


- holo 


. £l2£34---£56£78 

'£l3£l4---£57£58 


)l/4 


~ e, 


(3.31) 
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in the limit e —>• 0. Totally, there are four terms that are of order 0(e 4 ) in (3.29). 
Using (3.6), we get Afor n = 2, 

? (2) ^ r^(?l 2 ^34?56?78_)^^l/8 ^l)U))uHP 


A Sy ~ - log ( 


= - log ( 


(2e) 8 

. (^12^34^56^78)^ 


((02(6,6)02(6, 6))c/hp) 2J 


)l/8. 


(46 2 ) 1 / 2 


= 0, 


(3.32) 


(2e) 8 ' (£i2£34£ 56 £78) 1 / 4 J 

In L — l < t < L + l with e -+ 0, (3.12) shows £14 ~ £23 ~ e. Terms making non-vanishing 
contribution to 4-point correlation function will change in the limit e -+ 0. For example 
the third term in (3.29) as a sub-leading term is, 

{ 0 \ (£1, £i)Oi (£2,6)Oi (£3, £3)0! (£4, £4 ))uhp 

= (C ) l(£l ) £l)Ol(£2,£2)Ol(£3,£3)Oj(£4,£4)C ) i(£5,£5)C ) i(£6,£6)C ) l(£7,£7)0|(£8,£8))ij2 

,£i4£23---£58£67n 


- holo 


= ( 


-) 1/4 ~0(1), 


(3.33) 


'£l2£34-"£56£78' 

where ‘...’ stands for the terms that are 0(1) in the limit e -+ 0. One could count the 
leading contribution term by term in (3.29), there are only two terms that are of 0(e~ 1 ). 
Thus the variation of the Renyi entropy AS ' 1 ;”' 1 for n = 2 is 


A ~ — log 


= - log 


(2e) 8 

(£l4£32£56£78) 2 


)l/8 

)l/8. 


(£l4£32£56£78) 2 n 1/8 (^2 (£l> £l) ■ ■ 0 2 (£4, £4)) UHP 


((02(£i,£i)0 2 (£2,£2));7hp) 2 - 1 


(4c 2 ) 1 / 2 


= log 2 , 


(3.34) 


(2e) 8 ' 2 (£i 4 £ 2 3 £ 56 £ 78 ) 1 / 4 - 

2. For the Dirichlet boundary condition, i.e., 77 = —1. 

All the terms in (3.29) are non-vanishing in term of the neutrality condition. However, 
there are still four different terms in (3.29) contributing to the leading order in the region 
L + l < t ot t < L — Z in e —>• 0 limit. Thus the ratio (3.6) is 1 and A = 0. In 
the region L — l < t < L + l, there are two terms in (3.29) with e —>• 0 limit, which is 
the same as the situation of Neumann boundary condition. Thus the ratio (3.6) is ^ and 
A S y A = log2. In this example we could see that A S^' does not depend on the choice of 
boundary conditions. 

3.3 n-th Renyi Entropy for free boson 

In this subsection, we would like to generalize our studies to n-th Renyi entropy which 
involves in the the 2n-point correlation function on B n . The conformal transformation 
(3.4) (3.5) can map B n to UHP. Finally one can calculate the 2n-point correlation function 
on UHP by the ‘method of images’ in terms of 4n-point correlation function on R 2 . The 
points £ 1 , £2 ■ "£ 2 n on B n are 


£2/2+1 = ~i 


^ 2 mk/n +l ^l/n _|_ y 


£2/2+1 — —i 


^liirk/n j ^ 1 /n y ’ 

C^ink/n ^ 1 /n | y 

p 2 ink/n( £>1+^ \ 1 /n 1 ’ 

k uj\—1 / 


£2/2+2 = ~i 


^ink/n^W2+l ^1 In _|_ y 


g 2 ink/n _ y 


£ 2 / 2+2 — —* 


g 2 iirk/n'jl/n _|_ y 


p2i'Kk /n /In _ 1 ' 

' LdO —/ 2 


(3.35) 
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where 0 < k < n — 1. In the region t>L + lort<L — Z, one can obtain 

g ie 2iirk/n k 

S2fc+1 ~ &k+2 
&k+l — ?2fc+2 — 


n(uj\ - 0 1_1/n (wi + 0 1_1/n [e 2i7rfc/n (wi - 0 1/n - ( w i + l) 1/n } 2 

-8 ie 2ink / n le 


n 


{oj 1 - Z) 1-1/n (wi + /) 1 - 1 /n[e 2i7rfc / n (w 1 - Z) 1 /™ - (u)i + iy / n ]' 


and the Jacobi factor of conformal transformation, 

d&k+i d^k+2 1 t 

J - — J - — 7T~&k+l,2k+2 

d^2k+l d(jJ2k+2 2ze 

d&k+i d^2k+2 1 £ 

JZ - — 7 : - — ——?2fc+l,2fc+2 

dui2k+l dUJ2k+2 2ie 


(3.36) 

(3.37) 

(3.38) 


In L - l <7 t 7 L + Z, one could find 
£2 k ~ &k+i — 

? 2 fc+l - & 2 fc — 


8ie 2ink/ n i e 


n(—l + 0 Jiy~ l / n (l + wi) 1 - 1 /"[e 2i7rfe / n (a;i - Z) 1 /™ _ ( Wl + Z) 1 /™] 2 

-8 ie 2M / n le 


n(u>i - Z) 1 " 1 /™^ + _ Z)l/ n - (wi + Z) 1 /”] 

and the Jacobi factor of conformal transformation, 


d&k d&k+i 1 ^ 

j - — 7 -— Tr~^2k,2k+l 

dU2k duj2k +1 2ze 

d£,2k d^2k+l 1 7 

dco 2 k du 2 k +1 2ze 


(3.39) 

(3.40) 

(3.41) 


3.3.1 Local excitation 0\ 

We consider the operator 0 1 firstly. In the region L — l < t or t > L +1 with e —>• 0, the 
2n-point correlation function of Oi is 

(0}(6!?l)0l(6>6)---0l(6n,6n)) UHP 

= (0}(Cl,Cl)---0l(6n,6n)0j(^2n+l,6n+l)---0i(^ 4n ,^4„))ij 2 

— (0|(Cl,6)0l(?2 ; 6))---(0|(C4n-l,6tn-l)0i(^4 n ,^ 4 n ))ij 2 

(^12---C2n-l,2n) 1 / 4 (^2n+l,2n+2---C4n-l,4n) 1 / 4 ’ ^ ^ 

where 0\ is dehned as e - *^/ 2 with parity transformation and £ 271+1 = £1 v;£ 4 n = £ 2 n- 
Using (3.6), the variation of the n-th Renyi entropy can be obtained as follows. 

du: k h du k h ( 0 + (£1, £1) O (£2,£2) • • • 1 ( 6 n-1, ^ 2 n-l) 0 (^ 2 n, &n))uHP 


- holo 

. holo 




((Ot(wi, wi)0(w 2 , w 2 ))bi ) r 


1 — n 


log 


(Cl2---C2n-l,2n) 2 (C2n+l,2n+2---C4n-l,4n) 2 \ 1 / 8 


(2e 


1 4 n 


(4e 


2\n/4 


(^12---?2n-l,27i) 1,/4 (?2n+l,2n+2---C4n-l,4n) 1//4 - 


= 0 , 


(3.43) 
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where we have used (3.19)(3.36)(3.42), and Jacobi factor (3.37)(3.38). 

In the region L — l < t < L +1 with e —>• 0, the 2n-point correlation function on UHP 
is 

(0l(6,6)0l(6,6)---Ol(6n,6n)) UHP 
1 


Then, 


ASj 0 = 


(£23 • • • l ^2n,l) 1//4 (6n+l,2n+2 • • -^4n-l,4n ) 

' (?23-"?2n,l) 2 (?2n+l,2n+2---C4n-l,4n) 2 \ V 8 


(3.44) 


1 — n 


log 


(2e 


\4 n 




(4e : 


2W4 


(^23---^2n,l) 1 / 4 (^2n+l,2n+2---^4n-l,4n) 1 / 4 - 

where we have made use of (3.44) and Jacobi factor (3.13)(3.14). 


= 0, 


(3.45) 


3.3.2 Local excitation O 2 

For operator O 2 , there are 2 2ra terms making contribution to the correlation function like 
the ones in (3.29). Firstly let us consider the case with Neumman boundary condition. 
These terms with equal number of 0\ and 0\ can survive in the limit e —» 0. The 2n-point 
correlation function of 0\ on B n can be expressed by 4n-point correlation function on R 2 , 

(01 (6 5 6 ) Oi (£ 2 , 6 ) ■ • ■ ■ 0{ (£ 2 In-1,6 n- 1) Ol (6n, 6n )) UHP 


(0l(6,6)"-0l(6n,6n)Oi(6n+l,6n+l)---0l(6n,6n)) 


it 


- holo 


(3.46) 


To be convenient, we use the symbol +1 referring to ()\ and —1 referring to 0\ in the 
correlation function to simplify our analysis. Then the 2n-point correlation function on 
UHP can be formally written as 


(02(£l> £l)02(6,6)---02(6n-l,6n-l)02(6n,6n)){/HP 


- E 

on Z_^ 


(U) *2 j ■ ■ -U2n, *2n+1 ■ ■ ■ j 6n ) R 2 ■ 


(3.47) 


il,l2,...,*2n = ±l 

*2n + l =»1, • ■ .*47i=*2n 


where we have made use of the image method in the second line and ij = ±1 stands 
for the operator 0\ or 0\ with the coordinate (6,6)- The constraints i 2 n+j = ij with 
1 < j < 2n corresponds to the Neumann boundary condition in terms of image method. 
For the correlation function in the Neumann boundary condition case, the non-zero terms 
in (3.47) should satisfy neural condition 


*1 + *2 + + * 2 n + ■•■* 4 n — 0 . 

In the region t<L — lort>L + l with the limit e —> 
contribution of the 2n-point correlation function (3.47) are 

{i\ 16 ) 6 n; 6 n+1 •••, *4n)/J 2 


(3.48) 

0, due to (3.36), the leading 


E 


*l+*2 =0,43+14=0 
...i4n-l+Mn=0 
*2n+l=U...*47i=*2n 


— 2 n (ii, 12)^2 ...(*2n-l, *2n)i£2 (*2n+l, *2n+2)j? 2 ■ (*4n-l, Un)R? , 


(3.49) 


15 - 










there are 2 n terms that are leading divergence after considering the constraints, and these 
terms are all equal to each other. Equivalently, (3.47) can be written by the notation O\ 
and 0\ as 


2n (Ol(6)Ol(6))i? 2 ---(O{(6n-l)0l(6n))i?2...(O{(6n-l)0l(6n))i?2 

2 n 


it, 


it, 


(£l2 • ■ -6n-l,2n) 1//4 (6n+l,2n+2 • ■ • ? 4 n-l, 4 n) 1//4 

Using (3.47), the 2n-point correlation function on P> 2 n is 


( 3 . 50 ) 


(02(6)6)02(6,6)---02(?2n-l,6n-l)C ) 2(6n,6n)){/HP 

1 

(62 ■ ••6n-l,2n) O 4 (6n+l,2n+2- • • 6 n-l, 4 n) 1//4 


it. 


(3.51) 


which is same as the 2n-point correlation function of Oi in t < L — l or L + l < t. Thus 
the variation of the n-th Renyi entropy in the CFT with Neumman boundary condition is 


A S% ] = 0. 


(3.52) 


In L — l < t < L + l with the limit e —>• 0, due to (3.39) the leading contribution of (3.47) 
should satisfy following constraints 

6ra+l T 6ri+2 — 0 , 6ra+3 T 6 n +4 — 0 , ..., 6n—1 T 6n — 0. ( 3 . 53 ) 

Combining with im+j = ij and (3.53), one can obtain 

6+6 = 0 , 6 + 6 = 0,..., 6n-l+6rc = 0 . (3.54) 

In terms of (3.39), the leading terms of (3.47) should also satisfy following constraints 


6 + 6 = 0, 6 + 6 = 0,..., 6n + 6 = 0. 

With these constraints (3.53)(3.54)(3.55), the correlation function is 

('0^(6,6)02(6,6)-"02(6n-l,6n-l)02(6n,6n))[/ffP 

2 1 


( 3 . 55 ) 


2 n ( 63 -"6n,l)( 6 n+l, 2 n+ 2 • •• 6 n-l, 4 n) 1 ^ 4 


(3.56) 


Putting (3.56) into (3.6) with considering Jacobi factor (3.13)(3.14), the variation of the 
n-th Renyi entropy is 


A5? ) = log2. 


(3.57) 


In CFT with the Dirichlet boundary, the only difference with the Neumann boundary 
condition is the constraints 6 n+j = ij —> 6 n+j = —ij■ One could check that the leading 
order correlation function is same as the Neumann boundary condition. Thus the n-th 
Renyi entropy is not dependent on the choice of boundary condition. We do not repeat 
here. 
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3.4 Renyi Entropy in Ising model 

It is natural to ask how about the Ising model which is simplest unitary minimal model. 
There are three kinds of primary operators, i.e., the identity I, the spin operator <r and the 
energy operator e. There are also two kinds of parity transformation which involve in two 
kinds of boundary conditions, which correspond to two different parity transformations. 
One of the parity transformation [ 22 ] is 

a(z, z) —>• a(z, z), /j,(z, z) —>• /j,(z, z), (3.58) 

where p is the disorder operator. The other parity transformation [22] is 

a(z, z) —>• n(z, z), fj,(z, z) —>• a(z, z). (3.59) 


We would like to study the local excitation by the spin operator a with conformal dimension 
(h = jQ,h = j^j) in the same setup given in section (3.1). The variation of the Renyi entropy 
for the subsystem — l < x < 0 is given by (3.6). 

The two-point correlation function 


{a\u 1 ,u 1 )a{uj2,u;2))B 1 = r) (£i> CiM^ 

i =i d ^i d Zi 

f[( (3.60) 

Li d Zi dZi 


where £3 = £ 3*,£4 = £v,*,£] = iuii, if is the held with parity transformation. 

The 2-point correlation functions on UHP have already obtained in literature, e.g., 
[20] [23] [24], 




/ £l3^24 

^12^23^14^34 


M(*'), 


/ Cl3^24 

^12^23^14^34 




with conformal blocks [25] 


F + (x) = \J\J 1 — 3 / + 1 + \J Vl — x’ — 1 , and 
F _(.x ) = \j \/l — x 7 + 1 — \J y/\ — 


(3.61) 


(3.62) 


where x is the conformal cross ratio x = ^ 12 ^ 34 /^ 13 ^ 24 - F+( x ') an d E_(a/) correspond to 
different boundary conditions respectively. The leading behavior of the 2-point correlation 
function in e —> 0 is 


(a\uJi,u 1 )a(uj2,uj2))B 1 ^ 


V2 

(4e 2 )V 8 ' 


(3.63) for both boundary conditions. 


(3.63) 
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In t > L + l or t < L — l with e —>• 0, the leading behavior of 4-point correlation function 


is 


(^(£1, £lM£ 2 , 6 )°^ (£3, £3) O' (£4, £4 ))uHP — (^(£1, £l)cr(£ 2 , &))uHp{&H l i 3 , £sM£ 4 , U))uHP 

2 

“ (£l2£56) 1 / 8 (£34£ 7 8) 1 / 8 ’ 


(3.64) 


where we have used the 2-point function of Ising model on UHP (3.61). In terms of (3.6), 
Ais 



-log 

-log 


h (£l 2 £ 34 £ 56 £ 78 ) 2 Nl/l 6 (° 2 (£l» £l) •• - a 2 (£4, U))uHP ' 
( 2e ) 4 ((02(£l,£l)02(£2,£2));7Hp) 2 - 

h (£l2£34£56£78) 2 Nl/16 (4e 2 ) 1/8 1 = 

^ (2e) 4 j (£i 2 £34£56£ 7 8) 1 / 8 J ’ 


(3.65) 


In L — l < t < L + l with the limit e —>• 0, the 4-point correlation function (3.64) 
on UHP can not be factorized directly. We also use the image method, which states 
that the 4-point correlation function on UHP can be expressed as linear combination of 
the holomorphic part of conformal blocks of the 8-point correlation function on the full 
complex plane, with coordinates £i...£s- As we know in 2 dimension full complex plane, 
there are (n — 3) independent cross ratios for n-point correlation function. In our case, 
there are 5 independent cross ratios. Thus the 4-point correlation function on UHP can be 
expressed by conformal blocks [25] 


(c rt (£li£l)c r (£2,£2)o' t (£3,£3)cj(£4,£4 ))uHP 
= ( t )^^2A b C b J 7 [b;x,x 1 ,X2,...,x 4 ], 

£l2£23£l4£34 , 


(3.66) 


where we define the 5 independent conformal ratios x, x\, ..., X4, i.e. x = £i2£34/£i3£24> x i = 

£l4£l2/£li£42; :c 2 = £ll£ 23 /£l 2 £l 3 > ^3 = £ 2 !£2i/£22£l4> x i = £32£3i/£33£24) b ill Conformal 
blocks T runs over all the intermediate conformal families and the coefficients A p are 
determined by boundary conditions. Here we define £j = £j. The conformal blocks [25] 
satisfy the fusion transformation under the braiding operation [26][27][28], i.e., 


J 7 [b;x,xi,x 2 , ...,x 4 ] = F bc F[c ; 1 - x, x 1: x 2 , x 3 , x 4 ], 


(3.67) 


F bc is the fusion matrix. Making the fusion is equal to £2 •H* £ 4 . In summary, the leading 
divergence of the correlation function is related to c = 0 as identity as follows 

(3.66) = V F hc ( ) l sA b C b F[c-, 1 - x, Xl ,x 2 , ...,x 4 ]. (3.68) 

, £l4£43£l2£23 

o,c 

In terms of (3.12), one can find leading contribution of express (3.66) as following form in 
L — l < t < L + l with e —>• 0 

(3.66) ~ i ? oo(fJ 1 ’(£i,£i)fT(£4,£2))c/pp(cT' l ’(£3,£3)cr(£2,£4))pi7P, (3.69) 
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where 0 stands for the identity operator. Then the variation of the Renyi entropy for n = 2 

A4 2) = - log Too = log y/2, (3.70) 

where we have used the fact Too = [29] in two-dimensional Ising model. Note that the 

Renyi entropy does not depend on the choice of boundary condition in the Ising model, 
since the leading behavior of F + is the same as one of Tk with e —» 0. 

One alternative way to understand the phenomenon is to make use of diagram repre¬ 
sentation of conformal block as fig. [3]. Since the behavior of coordinates does not change 



Figure 3. Fusion transformation in Ising model. Where i and i presented in the figure stands for 
the position Z{ and zt respectively. Too[u] denotes the fusion transformation and 0 corresponds to 
the identity operator. 


when L — l<t<L + l, we only need one time fusion transformation, which is different 
with [4] . In the fig. [3] , 


Too [cl = T} 


oo 


a a 
a a 


(3.71) 


In L - l < t < L + l with e —>• 0, the leading contribution of (3.66) originates from the 
above conformal block involving in identity operator. Because the fusion factor Too can 
not be canceled in the ratio (3.6), A= — log Too = log \/2 [29], 


3.5 Renyi Entropy in General Rational CFTs 

In this subsection, we would like to generalize the analysis to the rational CFTs with a 
boundary in our previous set-up shown in subsection (3.1). In terms of (3.6), we should 
know the 2-point correlation on B\ and 2n-point correlation fuction on UHP as usual. In 
generic rational CFTs with a boundary, the 2-point correlation can be expressed as 


(O^wi, wi)0(cj2, w 2 ))bi = (0\Ci^i)0(^,2^2))uhp 

= ( 3 - 72 ) 


where A b are constants which are determined by the boundary condition, £ 3 , £ 4 , £ are 
given in (3.19). 4-point correlation function [23][24] of a primary function O on the R 2 , 
which can be expressed by 


(0\z 1 , z 1 )0(z 2 , z 2 )0 Jf (z 3 , z 3 )0(z 4 , z 4 ))r 2 = ^ 

b 


1 

(zi 3 Z24) 2h 


C b R[b ; 2 ] x c.c. 


(3.73) 
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In the limit z —>• 0 


T[b-,z\= z hb ~ 2h + ..., (3.74) 

where stands for higher order terms and z = Zl2Z34 and hh is conformal dimension of 

0 213^24 

primary class b. 

In terms of (3.74), (3.72) with taking £ —> 0 is 

1 —oh 

(0*(ui,u:i)0(w 2 ,w 2 )) Bl = (3.75) 

VS13S24I ) 

In t > L + l or t < L — l with e —>• 0, the 2n-point correlation function on UHP is 

(<^( 6 , 6)0(6,6)---O t (6n-l,6n-l)0(6n,6n))[//fP 
— {0\^i,E,i)0(^2,^2))uHP---{0\^2n-l,i2n-l)0{^2n,i2n))uHP- (3.76) 


The 2-point correlation on UHP is 

(o»« 2k-1^2k-l)0(^2k,^2k))uHP — 


A°C° 


(£ \2k-l,2k^2n+2k-l,2n+2k) 2h ’ 


(3.77) 


where 1 < k < n, 6 n+ 2 fc-i = 1- 

Taking (3.76)(3.77) (3.36)(3.72)and Jacobi factor (3.37)(3.38)into (3.6), one can obtain 


AS^ = 0 . 


(3.78) 


In L — l < t < l + L with e —>• 0, the 2n-point correlation function on UHP could 
be written as a linear combination of the holomorphic conformal blocks of the 4n-point 
correlation function. We make the following fusions [26] [27] [28] 


(£l > £ 2 ) (^3 1 £ 4 ) ■■■ (^4n—1 > ^4n) (?2,£ 3 ) (6, £ 4 ) ■ (?4n-l, ^4n) 

(6,6)(C4,?5)(?1>?6)--- ->• -t (^2,^3) — (Cl>C2n) —(C4n-l,^4n), (3.79) 

where £ 2 n+i = 0 < i < 2n. With using n — 1 times fusion transformation as the one 

shown in fig. [3], we get 

{O^ (Cl) 6)0(6; 6)-"O t (6n-l,6n-l)0(6n,6n))f/HP 

— ^ 00 [O] n-1 (Ot(£i, fl)0(^2n, &))uHP-~ (0^(^2n-2,^2n-l)0(?2n-l,C2n)) C/i/P, (3.80) 


We could calculate ASJ by using (3.6). The variation of the n-Renyi entropy is 

A = - log F 00 = log d 0 , (3.81) 

where do is the quantum dimension [29] of operator O. (3.78) and (3.81) are same as 
the case in rational CFTs living on the full complex plane. We apply the fusion rule of 
conformal block of 2n-point function on UHP to obtain the A S^K As we know, the 2 n- 
point correlation function on UHP can be expressed by linear combination of chiral part 
of conformal block associated with 4n-point correlation on full plane. In this subsection, 
we do not make use of parity transformation containing the boundary information. The 
boundary data has been encoded in the coefficient of conformal block in this subsection, 
i.e., A p in (3.73). We could see that AS\ also does not depend on the choice of boundary. 
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4 Conclusion and Discussion 


In this paper, we have investigated two kinds of effects on the locally excited states with 
time evolution. In the first case, we study the locally excited states with thermal effect 
in low temperature system. We have figured out the thermal correlation which is the 
same as [7] for the short interval limit. The Renyi entropy is equal to summation over the 
logarithmic of quantum dimension and thermal entropy in low temperature during the time 
L — l<t<L + l. In this paper, we just only confirm that such kind of sum rule is only true 
for the short interval l in the low temperature limit. We make use of different approach 
[6] to obtain the thermal correction to Renyi entropy which can be reduced to [17] in low 
temperature. One can also calculate the Renyi entropy in the large interval[7] limit, the 
higher temperature limit [30] as well as beyond the leading order of the perturbation (2.4). 
We expected the sum rule relation is still hold in those cases. But one should note that 
we actually do not consider the back-reaction of the locally excited states to the thermal 
environment. When the energy of the local excitation is much lower than the thermal 
environment, it is safe to ignore the back-reaction. But in some special situation we expect 
the sum rule will break down. It is an interesting topic to consider in the future. 

In the second case, we have studied the Renyi entropy of local excited states in 2 
dimensional CFTs with a boundary. For 2D CFTs with a boundary, to obtain Renyi 
entropy can be converted to obtain the correlation function on UHP by using of conformal 
transformation technique. As a warm up, the Renyi entropy has been calculated with 
help of image method in the 2D free field theory with a boundary. The Renyi entropy is 
vanishing for operator O 2 (3.15) in t>L + l or t<L — l and log 2 in L — l < t < l + L, 
which is the same as previous study [4] in full complex plane without boundary. To confirm 
this fact, the Renyi entropy have been calculated in Ising model and more generic rational 
CFTs. Although the correlation function and conformal blocks in 2D CFTs with boundary 
are totally different from the ones in 2D CFT without boundary, we get a same maximal 
value of the Renyi entropy for the rational CFTs without a boundary[4]. In generic 2D 
rational CFTs with a boundary, we confirm that the maximal value of Renyi entropy is the 
same as the one in 2D rational CFTs without boundary. [31] also try to understand the fact 
which is not contract with that the left- and right-moving chiral sectors are decoupled. [31] 
generalize the result in [4] to irrational CFT, for example, Liouvile CFT. They found that 
the left-right entanglement entropy saturates the Cardy entropy. In terms of standard view, 
the Cardy entropy counts the microscopic entropy of actual CFT spectrum. The Cardy 
entropy seems to suggest two chiral sectors are decoupled. The authors in [31] proposed a 
pragmatic point of view to reconcile [4] with the fact that there should also be a comparably 
large EE between the two chiral sectors of CFT. For example, Non-chiral local operators 
will be left-right entangled. In BCFTs, the two chiral sectors are no longer independent. 
We have shown some additional examples to confirm the pragmatic perspective. 

For general rational CFTs in 2D, the Renyi entropy highly relys on the conformal 
blocks of the theory. The n-point correlation functions in 2D CFTs with boundary are 
related to the holomorphic part of conformal blocks of the 2n-point correlation functions 
on the 2D full complex plane. This relationship had been studied by the image method 
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[20] [21] very well. More precise relation is that an n-point function in the UHP, which is 
a function of the coordinates (zi, ,z n ;z\, ...,z n ) behaves under conformal transformations 
in the same way as the holomorphic factor of a 2n-point function in the full plane which 
depends on (zi ,..., z n ; z*,z*), analytically continued to z* = Zj. In [4], the time evolution 
of Renyi entropy highly depends on the holomorphic part of conformal block. In 2D CFTs 
with a boundary, the boundary changes the evolution of the Renyi entropy but does not 
change the value of the Renyi entropy, which is closely related to fusion constants in the 
bulk. Because the behavior of the ‘image’ coordinates (anti-holomorphic coordinates) does 
not change as the holomorphic coordinates when L — l < t < L + 1, we get the same Renyi 
entropy as [4]. 

The boundary introduced here works as the infinite potential barrier for the time 
evolution of the entangled quasi-particles pairs [32] [3] [4] triggered by local excitation as 
shown in fig. [4]. The the Renyi entropy measures the entanglement between the quasi- 


0<t<L-l 

1 L-l<t<L 

1 L<t<L+l 

J 






X -L -1 A 

* -L -1 A 

X -L -/ A 


Figure 4. To show the time evolution of the entangled quasi-particles pairs. The boundary x = 0 
works as reflector. These red points stand for quasi-particles and the blue wave lines correspond to 
the entanglement between them. 

particles generated by local excitation. After entangled pairs are created at — L, the two 
quasi-particles will propagate in two opposite directions, i.e., left-moving and right moving. 
When the right-moving particle enter the interval —l<x< 0 denoted by A, the Renyi 
entropy takes maximal value due to entanglement between two entangled particles. In 
fig. [4], the blue wave lines mimic the entanglement of two entangled quasi-particles. When 
the right-moving quasi-particle reaches the boundary, the quasi-particle will be reflected 
by the boundary without losing energy. As the calculation in [33] shows the locally excited 
states carry the energy of 0(e~ 1 ). The conformal transformation, z —>■ z + e(z ) and z —» 
z-he(z), should keep the boundary conformal invariant, which lead to the constraint T = T 
on the boundary. In the Cartesian coordinates, the constraint becomes T x t = 0, which 
means that no energy can flow across the boundary. This is main reason the quasi-particle 
must be reflected no matter what is the conformal invariant boundary condition. In this 
sense, the boundary change the time evolution of Renyi entropy.In our paper we take the 
scale e as the minimal scale, and keep the leading order of e in the calculation. So we must 
miss some information when t ~ L, i.e., the quasi-particle is close to the boundary. In this 
case, we should make use of bulk boundary correlation functions and boundary structure 


- 22 - 








constants in BCFT to figure out the time evolution of entangled quasi-particles. The next 
leading order calculation of e may give us more insight on this point. 
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